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Abstract: We use the matrix model — gauge theory correspondence of Dijkgraaf and 
Vafa in order to construct the geometry encoding the exact gaugino condensate superpo- 
tential for the J\f = 1 U (N) gauge theory with adjoint and symmetric or anti-symmetric 
matter, broken by a tree level superpotential to a product subgroup involving U(Ni) 
and SO(Ni) or Sp(^) factors. The relevant geometry is encoded by a non-hyperelliptic 
Riemann surface, which we extract from the exact loop equations. We also show that 
0(1/N) corrections can be extracted from a logarithmic deformation of this surface. 
The loop equations contain explicitly subleading terms of order 1/N, which encode 
information of string theory on an orientifolded local quiver geometry. 
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1. Introduction 



Holomorphic terms in the effective action of supersymmetric gauge theories give much 
desired exact information about their non-perturbative dynamics. Based on geometric 
engineering through string compactifications, Dijkgraaf and Vafa [1, 2, 3] conjectured 
that the exact holomorphic superpotential of certain effective M = 1 gauge theories 1 is 
calculated by the planar diagrams of a matrix model. In the simplest case, an M = 2 
U(N) gauge theory is broken by a tree-level superpotential W((f>) = X^=i ^f~j"Tr0 fe+1 
to an Af = 1 theory with gauge group G = Y\i=i U(Ni). Here is an adjoint chiral 
superfield, and one assumes that all zeros of W are simple and that YH=i N% = N. 
An effective theory emerges by integrating out the massive chiral fields. Its effective 
low energy Lagrangian can be expressed in terms of the gaugino bilinear fields Si = 
— ■^2TrW ai W ai where % indicates the factors of the unbroken gauge group. The claim 
is that the effective action as a function of these fields is calculated exactly by the 
leading terms in the 1/N expansion of a matrix model whose matrix potential is given 
by W. 

Technically the field theory results for the case mentioned above are not new. They 
were obtained by geometric engineering in [5] and partly by factorization of the M = 2 
Seiberg-Witten curve [5, 6]. This allows for explicit checks [7] and leads to a particular 
interpretation of the matrix model. The relevance of matrix model structures, such as 
the planar loop equations, can also be understood more conceptually using supergraph 
techniques [8] or anomalies [9, 10]. Further developments can be found in [11, 12, 13, 14]. 

In this paper, we shall add matter in the symmetric or antisymmetric representa- 
tion of the U(N) gauge group and consider an additional term Tr Q<pQ in the tree-level 
superpotential 2 . This leads to a richer vacuum structure, which has branches where 
additional SO(Nj) and Sp{N i /2) gauge factors appear. The conjecture of [1] then gives 
us a rough recipe for how to obtain a dual matrix model. However, supersymmetric 
vacua are described by a symplectic quotient, which is equivalent with a holomorphic 
quotient of the space of F-flat configurations through the action of the complexified 
gauge group. Thus the gauge-theory analysis of the vacuum structure leads to solutions 
in the space of complex matrices modulo complexified gauge transformations. On the 
other hand, Hermitian matrix models have hermiticity constraints on their matrix vari- 
ables. To make the matrix model useful for the gauge theory analysis, an appropriate 
constraint should be imposed in a way which allows one to recover the correct vacuum 

1 Some of these results carry over to N = orbifolds [4] of N = 1 gauge theories. 

2 Symmetric and antisymmetric representations are interesting because they are building blocks 
of chiral supersymmetric gauge theories. In particular many models with dynamical supcrsymmetry 
breaking are based on antisymmetric representations, see e.g. [15, 16, 17]. 
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structure of the gauge theory. As an example, consider the one-matrix model associated 
with the U(N) gauge theory broken by a cubic superpotential W((f>). Then the field 
theory has two classical vacua at W'{4>) = 0, which, for a Hermitian matrix 0, would 
correspond to a local maximum and a local minimum for its eigenvalues. Perturbation 
theory around the local maximum would make no sense in the Hermitian model. In 
this case, the relevant information was extracted in [18] and [19], where it was shown 
that it matches the gauge theory instanton expansion relevant for both vacua. Using 
the BRST formalism, the residual gauge symmetry U(N)/(U(Ni) x f/(AT 2 )) acting on 
the vacuum leads to an expansion scheme in terms of a two-matrix model interacting 
through the ghosts [7]. To avoid convergence problems, one must reinterpret the ex- 
pansion around the extrema in terms of a Hermitian N± x N± matrix interacting with 
an anti- Hermitian N 2 x N 2 matrix [18]. The resulting perturbative expansion then 
recovers the leading F-term gravitational corrections to the gauge theory, which in the 
matrix model arise at order 0(1/ 'N 2 ) [18, 19]. A more systematic analysis of such issues 
requires a holomorphic construction of the relevant matrix model [20] . Namely it leads 
to a so-called 'holomorphic matrix model', whose relevance was already pointed out 
in [1]. As shown in [20], the procedure employed in [18] is indeed perfectly justified 
in the holomorphic matrix model. In the present paper, we shall encounter similar 
problems even before choosing a vacuum. Accordingly, we must carefully choose a 
'real section' through the space of complex matrixes in order to make all integrals con- 
verge. Moreover, we shall need a complex regularization of the matrix model in order 
to avoid introducing spurious constraints on its filling fractions. This is very similar 
to what happens in the case of holomorphic A 2 models [20], where use of a complex 
regularization is crucial in order to avoid similar problems. 

In section two we introduce the gauge theory model and analyze its classical vacuum 
structure. Using threshold matching, we also derive the leading logarithmic terms in 
the effective superpotential. 

In section three we define the corresponding matrix model, whose action can be 
taken to be the superpotential of the gauge theory. As mentioned above, it will be 
crucial to choose a suitable replacement for the hermiticity constraint. Contrary to the 
one-matrix model, we cannot impose hermiticity of the matrix M which corresponds 
to the adjoint chiral multiplet <fi. Indeed, this requirement would prevent the matrix 
model from probing the complete vacuum structure of the dual gauge theory. As 
explained in [20], this is not a particular property of our model, but can be traced back 
to similar issues encountered in ADE matrix models. Using a holomorphic definition 
and a 'complex' regularization of our model, we shall extract two exact loop equations 
(Ward identities) which are cubic and quadratic in the resolvent. A new feature of 
these relations is that they contain explicitly terms of order 0(1/N) and 0(1/N 2 ). 
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Taking the large N limit of these constraints allows us to extract a non-hyperelliptic 
Riemann surface which governs the planar solutions of the model, as well as the effective 
superpotential of the dual field theory. The dual string geometry is expected to be 
given by an orientifold of the non-compact Calabi-Yau background dual to the A 2 
quiver theory 3 . We therefore also give the explicit derivation of the full loop equations 
for the matrix model based on the A 2 quiver, and relate it to the Z 2 orbifold and 
orientifolds thereof. However, there are subtleties in the orientifolding, in the large N 
transition and in the precise definition of the 5-model cycles. We therefore prefer to 
use matrix model techniques in order to extract the Riemann surface, which should 
already encode all relevant information about the holomorphic terms in the M = 1 
effective action. This surface is a triple cover of the complex plane, which cannot 
generally be written as a hyperelliptic curve. Since we have little guidance from a 
string geometry, we must intrinsically understand this surface, which in the geometric 
engineering approach would be obtained by integrating out two directions of a dual 
Calabi-Yau geometry [5]. In particular, one must take into account issues of non- 
compactness and relations between periods, which are important for a proper count of 
this curve's parameters. 

Similar to the case of SO(N) and Sp(f ) groups [21, 22, 23, 24, 25, 26, 27, 28, 29], 
the 't Hooft expansion leads to unoriented ribbon graphs and the large N expansion of 
the matrix model is spaced by N and not iV 2 . Thus the explicit 1/N terms entering 
our loop equations have field-theoretic relevance. To properly analyze such effects, one 
must implement the filling fraction constraints iVj =const of [1] in the finite N model, 
which will be done in Section 4 by introducing chemical potentials and performing 
a Legendre transform to extract a microcanonical generating function. This gives a 
direct derivation of a set of special geometry relations, and provides their finite N 
generalization. It also naturally leads to a finite N version of Whitham-type constraints. 

In Section five we show how 0(1/N) contributions can be computed from a mod- 
ified Riemann surface, which is obtained by performing a logarithmic deformation of 
the matrix model. We show that the 0(1/N) term of the microcanonical generating 
function can be computed by differentiating the leading (N = oo) contribution with 
respect to the coupling constant of such a logarithmic deformation. 

In Section 6, we check agreement between the matrix model and the field theory. 
For this, we use the matrix model to compute the leading (Veneziano-Yankielowicz) 
contribution to the effective superpotential for different branches of the moduli space. 
This is done by performing BRST gauge-fixing in the manner of [7, 18] and integrating 

3 Factorization of the M = 2 curve cannot give hints, since it leads to a completely different branch 
of the moduli space. 
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out the quadratic terms in the action. The comparison to the O(N) model and the 
A 2 -quiver is summarized in the Appendix. 



2. A first view of field theory properties 

In this section we study U(N) M = 1 supersymmetric gauge theory with matter in 
the symmetric or in the antisymmetric representation. The starting point is an Af = 2 
U(N) gauge theory with matter in the symmetric or antisymmetric representation 
(Q T ,Q T ) = s(Q,Q). We choose a tree-level superpotential 



which softly breaks the Af = 2 supersymmetry to Af = 1. In many of our considerations 
below it helps to view the theory as coming from a Z 2 -orientifold of the A 2 quiver 
theory. Orientifolds of supersymmetric A 2 quiver gauge theories have been constructed 
before with the help of Hanany-Witten type brane configurations in Type IIA string 
theory [30, 31]. These models were further analyzed in [32, 33, 34, 35]. Whereas an A 2 
quiver gauge theory contains two independent unitary gauge groups and matter that 
transforms under the bifundamental representation, the orientifolded model identifies 
the two gauge groups. In addition the orientifold modifies how the matter content 
transforms. In the A 2 model the chiral superfield Q transforms as U\QU\ where XJ\ 
and U 2 are the gauge transformation of the two independent gauge groups. In the 
orientifolded model this changes to Q — > U QU T . Therefore Q transforms in a two-index 
tensor representation of U(N) and the two possibilities of symmetric or antisymmetric 
representation correspond to the two choices of orientifolds. 

In the following we describe the vacuum and its phase structure of the emerging 
Af = 1 theory and calculate in various phases the Veneziano-Yankielowicz potential 
from threshold matching. 

2.1 The classical moduli space 

We are going to classify the possible constant solutions of the classical field equations. 
Throughout the analysis we also assume that Af = 1 Fayet-Iliopoulos terms are absent. 
As is well known in this case the space of solutions of the field equations can be obtained 
by minimizing the superpotential W and dividing by the complexified gauge group. The 
field equations are 



W = tr {W{(f>) + Q<PQ) , 



(2.1) 



QQ + W\<t>) = , 
cj)Q + Qcj) T = , 
Q4> + <f) T Q = . 



(2.2) 
(2.3) 
(2.4) 
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For quiver gauge theories similar field equations arise and the classical moduli space 
has been analyzed in [36]. We will adapt the methods there to the case at hand. 
Note first that from (2.3), (2.4) it follows that 

0"Q = Q(-0 T y\ 

Q^ = {-fYQ, (2.5) 

We multiply (2.2) from the left with Q, from the right with Q and commute Q through 
W with the help of the previous relation. This gives QQ[QQ + W(— 4> T )] = 0. Now 
we transpose (2.2) and use it to eliminate QQ and we obtain — W'(4> T )[— W'(4> T ) + 
W'(— 4> T )} = 0. Finally we transpose this last equation and arrive at 

[W'{4>) - W'(-(f))}W'(<f)) = (2.6) 

Equation (2.6) will be solved for <fi being diagonal <fi = diag(ai.ljvi, • • • , a n .l Nn ) and 
Y^i=i Ni = N. The entries a, L have to be the roots of one of the two equations 

W'(x) -W'(-x) = 0, (2.7) 
W'(x) = . (2.8) 

The vev of breaks the gauge group according to U(N) — > Yli=i U(Ni). The field Q 
decomposes into Qij with Qij transforming as a bifundamental (JVj, Nj) under U(Ni) ® 
U(Nj) if % j and as symmetric (antisymmetric) if % = j. An analogous statement 
holds for Q. From (2.3) it follows that 

(a, + aj )Q t3 = (2.9) 

and thus Qij = unless dj = — Oj and Ni = Nj. Such pairs of solutions are indeed 
generated by the roots of (2.7). Another special root of (2.7) is x = 0. Let us study 
now in more detail the different solutions. 

2.1.1 Pair of solutions a { = -aj = b ^ with W'(b) = W'(-b) 
In the relevant subspace Q and Q have to be of the form 

A gauge transformation acts on Q as 



U x q\Jl 
sU 2 q T Ul 



fl-L^r/r ) (2-11) 
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Since U\ and U2 are independent GL(Ni,C) matrices we can bring Q into the form 



9=(,°J'l- (2-12) 



From (2.2) it follows then that 



Q = - w '( b )(° 1 S o) ■ ( 2 - 13 ) 

The gauge transformations that leave these matrices invariant are given by (Ui)^ 1 = 
and therefore the residual gauge group is U(Ni) in this branch. 

2.1.2 Solutions with a; = 

This is the special solution of (2.7). The unbroken gauge group in this branch is U(Ni) 
and it acts on Q as 

Q -> UQU T . (2.14) 

Let us consider first the case of symmetric Q. Since U is a general linear matrix we 
can choose U such that Q becomes the iVj-dimensional unit matrix Q = 1. It follows 
then that Q = — W'(0)1. The gauge transformations that are left over have to fulfill 
U.U T = 1 and are therefore elements of SO(Ni). 

In the case where Q is antisymmetric similar arguments show that Q can be brought 
into the form 

/e0...\ 



Q = 



1 
-1 



(2.15) 



e ... 

Q is determined by (2.6) and the unbroken gauge group is Sp(Ni/2). 
2.1.3 Solutions with W'(ai) = 

In this case the equations of motion imply Q = Q = and the unbroken gauge group 
is U(Ni) in this branch. 

2.2 The Veneziano-Yankielowicz potential 

After having established the structure of the moduli space of vacua we derive now 
the leading terms in the low energy effective superpotential. Let us take a generic 
vacuum with = diag(0jv o , ail Nl , . . . , a n l Nn , M^, -M^, • • • ,h^N k , -MjvJ- The 
unbroken gauge group is 



i=l j=l 
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At low energies the non-Abelian factors in the unbroken gauge group confine and leave a 
U (l) n+fc Abelian gauge group unbroken. The effective dynamics of the gaugino conden- 
sates in the non-Abelian factor groups is captured by the usual Veneziano-Yankielowicz 
superpotential. For the unitary group factors this is 

Wi Y = S j ]og[-p±j , (2.17) 
and for the orthogonal or symplectic factor group it is given by 

/ A 3(AT -2 S )\ 

"tr (218) 
Taking all gaugino condensates into account we have 

n+k 

W eff = J2 W vz + W vz, (2-19) 

where we substitute Nj = Nj_ n if j > n. 

The low energy scale can be determined by threshold matching at the scales of the 
masses of the various massive W^-bosons and matter fields. In particular we find for 
the vacuum with gauge group U (iVj) 

n k 

a£S,. = (v"^)) Nt - a >r 2N3 11^ - h2 iT™ l ^T 2No >< 

jfr 1=1 

x f[(a t + a,)* n(a? - b^ia^ia^A^ . (2.20) 

r+i t=l 

The first term on the rhs of this equation comes from the fluctuations of in the iVj'th 
diagonal block around the vev Oj. The next terms stem from the massive off-diagonal 
W bosons and in the second line we collected the contributions from the matter fields 
Q and Q. Analogously one can analyze the other scale matching relations. We find: 

A£?,i + n = (V"(b t ) + V"(- bl )f f[(% - a])'™* n(6? - b>)~^(2b t )-^(b t )- 4No x 

3=1 l^i 

n k 

x n<«? - b2 ) Nr - ^w^^r^Z'^ ■ ( 2 - 21 ) 

r=l t+i 
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and 



It ru It, ft, 

aIs _2s) = (^^o))^- s nKo-^n^)" 4 ^n^)" r n^)" i < g i 1 2s - ^) 

j=l 1=1 r=l t=l 

In (2.21) the factor of 2 in the exponent of A high reflects the diagonal embedding of the 
U (iVj) gauge groups 4 . The higgsing by the vevs of Q which breaks U (Nj)<S>U (Nj) to this 
diagonal U(Nj) and the U(Nq) factor group to either orthogonal or symplectic groups 
produces also some massive fields. These come however always in the multiplicities 
of iV = 4 multiplets and therefore do not contribute to the threshold matching. To 
be specific, consider the breaking to SO(N ). The massive fields include W bosons 
and components of that lie in the coset U(N )/ SO(N ). They transform under the 
symmetric representation of SO (No). In addition Q and Q are also symmetric in this 
case. Thus we count (in M = 1 language) ^-(Nq + 1) vector multiplets and 3^-(N + l) 
chiral multiplets that receive masses at the scale set by the vev of Q. 

With this method we can of course obtain only an approximation to the exact low 
energy superpotential. But is it is expected to capture the correct logarithmic behavior 
at Si — 0. In section 6 we will derive this part of the superpotential from a one-loop 
calculation in a matrix model. 

3. The matrix model 

In this section, we construct the matrix model which is expected to calculate the exact 
superpotential of the M = 1 supersymmetric U(N) gauge theory with matter in the 
adjoint and symmetric or antisymmetric representation. Using direct manipulations, 
we shall derive the exact loop equations of this model. These are two independent 
identities, which are respectively quadratic and cubic in lv(z) and ou(—z), where ou(z) 
is the resolvent of the model. The large N limit of these relations gives algebraic 
constraints on the planar limit of the averaged resolvent, which lead to a proposal for 
the algebraic curve governing the dual gauge theory. We also give evidence that the 
planar vacuum structure of the matrix model agrees with the field theory. 

4 For the threshold matching of the SO/Sp factor group there arises the following well-known 
problem [37]. Projecting SU(N) to SO(N) or Sp(^) there are some roots which are invariant under 
both projections. However, these roots serve as long roots for SO(N) but as short roots for Sp(^). 
This results in a relative factor of two in the normalization of the roots in the projected groups and 
this normalization influences the indices of the representations, e.g. for the adjoint representation 
Cq = (9, 9)g y , where 9 is the highest root and g y the dual Coxeter number. 
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3.1 Construction of the model 

Our matrix model results by performing one of two orientifold projections on the matrix 
model proposed in [2] for the A 2 quiver field theory. The discussion of the covering A 2 
quiver theory and that of a related orbifold model can be found in Appendix A. 

We use the superpotential (2.1) of the M = 1 supersymmetric gauge theory with 
gauge group U (N), a chiral superfield (f> in the adjoint representation, a chiral superfield 
Q in the symmetric or antisymmetric two-tensor representation and a chiral superfield 
Q in the corresponding complex conjugate representation. In order to define the matrix 
model we have to choose a real section in matrix configuration space. In the one matrix 
model this is usually done by mapping $ to a Hermitian matrix M (as discussed in 
[20], this prescription is justified if the relevant superpotential is a polynomial of even 
degree). Following a similar prescription for our models would lead to a series of 
problems. Similar to what happens in the quiver matrix models [38, 39, 40, 41] and the 
O(N) [42] we would then have to require that all eigenvalues of M lie on the (strictly) 
positive real axis. However, we saw in the analysis of the classical moduli space that 
vacua with negative or zero eigenvalues play an important role. So it seems quite 
unnatural from the gauge theory point of view to restrict to such matrix configurations 5 . 
Instead, we shall consider the matrix model: 

Z N ,s,s = Vol ^ (jV) x J dMdQexp [— TVtr V(M) + itr Q^MQ] , (3.1) 

where the M-integration is performed not over the set of Hermitian matrices but rather 
over the set: 

M = {M e Matjv(C) | M- M* = 2iel} , (3.2) 

where e is a small positive quantity. We have also imposed the condition Q = —iQ^. 
Together with the shift of the eigenvalues of M into the upper half plane, this renders 
the Q-integration finite without further restricting the range of M's eigenvalues. Then 
the Q-integration is performed over the set: 

Q={QE Matjv(C) I Q T = sQ } , (3.3) 

where s = +1 and s = — 1 distinguishes the symmetric from the anti-symmetric repre- 
sentation. The measures dM, dQ are defined through: 

dM = n dMii n dRe dim M ^ ' 

i i<j 

dQ = Y[ d Re Qij d Im Q tj ( JJ d Re Q lt d Im Q^j * M . (3.4) 

i<j i 

5 A similar issue arises in quiver matrix models, as discussed in [20]. 
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This model has the gauge-invariance: 



M -> CMC 1 

Q -> f/Qf/ T , (3.5) 

where C/ is an arbitrary unitary matrix 6 . Note that we have included the inverse volume 
of the gauge group into the definition of the partition function (3.1). Also note that we 
are rather valiant about the convergence of (3.1), which is assured — with our choice 
of integration manifold — only provided that W is a polynomial of even degree. As in 
[20] , a consistent construction for odd degree potentials would require that we constrain 
the eigenvalues of M to lie on a certain path in the complex plane, whose asymptotic 
behavior is determined by the leading coefficient of W. In the present paper, we shall 
ignore this and related issues, a complete treatment of which requires the full machinery 
of [20]. 

Performing the Gaussian integral over Q yields 

fa**™ . H ~nr^(nir (3-6) 



where Ai, . . . , \n are the eigenvalues of the matrix M. 

Next we rewrite the integral over M in (3.1) in terms of eigenvalues of M: 



with 

N 



/ | J d\ k ' 

JR+ie k=1 



Cn,s Z N)Si£ , (3.7) 



N 2 S A n 2 N /2+N (■ \N(N+s)/2 , , 



Vo\(U(l) N x S N ) 

where the 'effective' action <S'(Ai, . . . , A at) in the eigenvalue representation is: 

k k k^^l k ,/ 



6 This can be seen as follows: The action and the measure c!M for the adjoint field are obviously 
gauge invariant. Therefore we concentrate on the part coming from dQ. It is easiest to switch for 
a moment to an index notation in which the fields transform as Qij — ► \{Ui Uj + slli 1 Uj k ) Qki, 

Qij Qkl 1 (Jjt i Jjj J + s jjt i Jjt We took 

care explicitly of the symmetry properties of Q and Q. 
The measure dQ picks up the product of the determinants of the transformation matrices defined in the 
previous equations. Since the product of the determinants is the determinant of the product we com- 
pute the product of the transformation matrices: \(U t k U 3 l + sU t 1 Uj k ).\{U ] k m U} n + sU} m U ] k n ) = 
2 (6™ SJ + s8f 8J 1 ) This is just the identity in the symmetric and antisymmetric representation respec- 
tively and of course the corresponding determinant is one. 
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In the limit e — > (which we shall ultimately take below), the factors e ln ( A fc+ A ') = 
Xk \ Xl in (3.8) produce well defined distributions according to the Sokhotsky formula: 

= V ( —i— )- i7i5(\ k + A,) , (3.10) 



A fc + A/ + i0+ V^fe + A 

where we took the limiting eigenvalues to be real in order to make the is prescription 
explicit and where V denotes the principal value. Substituting (3.10) into (3.8) (before 
exponentiating the Vandermonde factors) leads to an expression for the e — > + limit of 
Z as a sum over 'reduced' integrals. This gives a limiting statistical ensemble, which can 
viewed as the orientifold of the limiting ensemble extracted in [20] for the holomorphic 
A 2 model. 

The potentials V(x) we want to consider are of the form 

d t 

V(x) = t-i ln(x) + W (x) where W(x) = V -Ar x k+1 . (3.11) 

K~\~ -L 

k=0 

From (3.9) we see that the effect of integrating only over symmetric or anti-symmetric 
Q is a logarithmic correction to the potential of order 1/N. If we introduce 

U(x)=V(x) + ^\n(x) (3.12) 

the action S depends on s only implicitly through U(x). For the analysis of the matrix 
model it is useful to further introduce the following quantities. The expectation value 
of a function 0(\i, . . . , \ N ) is given by 

1 r N 

(°) = ^ ]]d\ k O(\ u ...,\ N )e- N2s (3.13) 

Z jR+ie k=1 

and the eigenvalue density p(A) and the resolvent cu(z) are defined to be 

pW = ^J2 5 ( x - x ^ and w W = ^Ejza^- (3 - 14) 

The two quantities are related by 

p(A) = hm "i^") ~ "V-i") , u{z) = f dX . (3 . 15) 



v^o 2%i ./r+k z-X 



13 



3.2 The microcanonical ensemble 

The framework of Dijkgraaf-Vafa requires that the model obeys certain filling fraction 
constraints. In [1], such conditions were imposed only on the large N microcanonical 
generating function, which is insufficient in our case since we will have to consider 
0(1/ N) corrections as can be seen from (3.12). This requires that we impose such 
constraints on the finite N microcanonical generating function, rather than on its large 
N counterpart. The relevant constraints are easiest to formulate by employing a mi- 
crocanonical ensemble. As we shall see below, the original path integral defines a 
(grand) canonical ensemble at zero chemical potentials. This allows one to recover 
the microcanonical generating function by introducing non-vanishing chemical poten- 
tials (which are canonically conjugate to the filling fractions) and then performing a 
Legendre transform to replace the former by the latter. 

3.2.1 The (grand) canonical partition function associated with a collection 
of intervals 

Let us cover the displaced real line R + ie with disjoint nonempty line segments A a 
subject to the condition: 

AiU---UA r = R + i£. (3.16) 

We shall let Xa denote the characteristic function of A a , and consider the filling fractions 
fa of A a : 

/« = ]7 • ( 3 - 17 ) 

k 

This gives the expression of f a in terms of uj(z): 

f a = J d\p(X)xa(X) = jf d\p(X) = j ^-uj(z) , (3.18) 

Picking chemical potentials fj, a , we consider the (grand) canonical ensemble associated 
to our collection of intervals 7 : 

„ N r 

Z(t,n)= ]JdX k e- N2s - where = S + ^ p a fa ■ (3.19) 

The original partition function results by setting /i a — 0. Introducing the (grand) 
canonical generating function: 

F(t,n) = -^\nZ(t,fi) , (3.20) 



7 In the following /U, S may stand for r-tuples and t for (t-i,to, ■ ■ ■ , td). 
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we have the standard relation 

d 



dfi c 

In this subsection, the brackets (...) always denote the expectation value taken in the 
(grand) canonical ensemble. Since the union of A a covers the whole integration range, 
the expectation values of the filling fractions fulfill the constraint: 

r 

£</«) = 1 . (3.22) 

a=l 

3.2.2 The microcanonical generating function 

Following standard statistical mechanics procedure, we define: 

S a ■= J-f (3.23) 
and perform a Legendre transform to extract the microcanonical generating function: 

r 

F(t, S) '■— ^ S a /i a (t, S) — F(t, fj,(t, S)) . (3.24) 

a=l 

In this relation, ji a are expressed in terms of t and S by solving equations (3.23). The 
constraint (3.22) shows that S a are related through: 



J2 S « = 1 > ( 3 - 25 ) 



a=l 

so we can take S 1 . . . S r -i to be the independent variables. Then equations (3.23) 
express /i a as functions of t and these coordinates, and equation (3.24) implies: 

dF 

fJ>a~ Hr= for a = 1, . . . , r - 1 . (3.26) 

oS a 

Note that /i a are only determined up to a common constant shift; this is due to the 
constraint (3.25) on S a . 

Working with F(t, S) amounts to fixing the expectation values of the filling fractions 
by imposing the quantum constraint (3.23): 

(/«> = jf ^-Mz)) = S a , (3.27) 

with S a treated as fixed parameters. This gives a meaning to the procedure of [1] 
beyond the large iV limit. 
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3.3 The quadratic and cubic loop equations 

The classical equations of motion following from the action (3.9) are dS/dXk = 0, 
where: 

N- 

' I \ i ' A 



I^W-iEA^Ei. (3-28) 



The action (3.9) describes a system of charged particles moving along M + ie, together 
with a set of mirror charges of opposite sign. For s——l, each particle interacts with 
all mirror charges excluding its own, while for s—1 it is also attracted to its own mirror 
image. 

The partial derivatives (3.28) will be used below to obtain a set of Ward identities 
for the matrix model. In turn, we shall use these identities in order to derive two loop 
equations for resolvent u(z). Consider the integral: 

w i L i A ? ^ ( MK "-^ e ~ N ' s ^ = (3 - 29) 

Differentiating shows that (3.29) is equivalent to the Ward identity 

^e<(^)4e^4e^)*-^> - ° 

3.3.1 The quadratic loop equation 

The quadratic loop equation is obtained by setting ip k = (z— A^) -1 . Together with the 
identity 

2 ^ ~^Xk + ^ {z-\ k ) 2 = ^ ~^X k (3 ' 31) 

one obtains the relation 

<»<*> -^E^-^E^^) = °- < 3 - 32 > 

Adding to (3.32) the same equation with z replaced by — z and absorbing the depen- 
dence on the potential into f(±z) leads to the quadratic loop equation: 

(co(z) 2 + uj(z)Lo(-z)+Lj(-z) 2 -U'(z)uj(z)-U'(-z)Lo(-z) + f(z) + f(-z)) = . (3.33) 

In deriving this we used the identity: 

1 1 1 1 V- 1 1 , nni . 
> h > > (3.34) 
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and introduced the quantity: 

/W .1 ? 5MZ^. ,3,5, 

In the analysis of the large N Riemann surface below we will need to know the properties 
of f(z) in more detail. From the definition of U,V,W in (3.11) and (3.12) it follows 
that f(z) can be written: 

f(z) = f(z)-<j>- t - 1 + s/{2N) where = 1 £ A," 1 = - w (0) . (3.36) 

k 

The function / (z) is a polynomial in z defined through: 

m __^^^> (33?) 

Since oo(z) ~ l/z as z — > oo, one finds that the large z behavior of / is /(z) ~ t d z d ~ x . 
In particular the polynomial f(z) has degree d— 1. 

3.3.2 The cubic loop equation 

The cubic relation for uj{z) results from substituting 



, so that: 



N ^ (X m +X k )(z-X k ) 

W k ~ N ^ (X m +X k )(z-X k ) \7^X~ k ~ Xm+Xj + N 4 X| (*-A fc ) ( ' } 
into the Ward identity (3.30). To proceed we make use of the relations: 



— T 



^ 3 (Afc-A,)(A m +A fc )(z-A fc ) 

m,k,l(^k) 



and 



z ) + iV 3E (Am+Afc)Mfc) C_ Afc A m +A fe ) 

fe,m 

+ A^^(A m +A fc )(A z +A fe )(^+A fe ) (3,39) 

AT ^ A| (*-A*) z> z N k ^k 1 J 
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Equations (3.39) and (3.40) follow from the partial fraction decompositions: 
1111 

+ 7 7T7 7T7 7T + 



(z— a) (z— b) (z— c) (b—a)(c—a)(z—a) (c—b)(a—b)(z—b) (a—c)(b—c)(z—c) 
1 1 1 1 , 

+ r, iSTI ^ ( 3 - 41 ) 



(2-a) 2 (2-6) (a-b) 2 (z-b) (a-b) 2 (z-a) (a—b)(z—a)' 

When substituting (3.38), (3.39), (3.40) into (3.30) some of the terms cancel and we 
are left with 

2 1 U'(X k ) 1 1 / 1 1 \ 

-'UJ w{-z) + Jf 2 l^ (Xm+ x k )( z -X k ) + Jp^. (X, + \ k )(X m+ X k ) ~ ^X k ) 

k,m 



(X m +X k )(z-X k ) N 3 , (Xi+X k )(X m +X k ) \z-X k -z 



k 



As in the quadratic case we can add to this the same equation with z replaced by —z. 
This removes the triple sum and gives 



uj(z) uj{ 



, X / N / N2 1 U'(X k ) ( 1 1 \ 



1 a;(z)+a;(-z)-2a;(0) \ _ 

iV^ 4^ / " U [6A6) 

To proceed we introduce the quantity 

As with /(z) in (3.36), for later use we note that g(z) can decomposed into a polynomial 
part g(z) and a pole 1/z as 

^(zj = — 7 • where 

1 W'(z) - W'(X k ) , 1 1 cj(0) 2 , n 

^) = ^E (A ,V J A fc ) ( A) md ^^^^ (3-45) 

As before, since W'(;z) has degree d the large z behavior of the polynomial g(z) is 
(const) • z^" 1 , so that g(z) has degree d—1. The expression 7 = u(0) 2 /2 has been 
obtained by evaluating the identity (3.34) at z—0. 

Using (3.32), one finds the following relation for the expectation value of g(z): 

<»W> = (-j^E ( ;JSU) + "*« (<"«' + / W " W*)) ) (3-46) 
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Substituting (3.46) into (3.43) we obtain the cubic loop equation: 

< u(z) 2 u;(-z) - g{z) + U\z) {u{zf + f(z) - U\z)u{z)) + (z <-> -z) ) 

1 (u{z)+u{-z)-2u>W) n ^ 
iV 2 Az 2 ' 1 J 

3.3.3 Loop equations in terms of contour integrals 

The two loop equations (3.33) and (3.47) can be presented in a more compact form 
when using contour integrals. This form can be obtained as follows. Let 7 be a contour 
that encircles all the eigenvalues but not the point z and not the poles of u>(—z). One 
can verify the two identities 

dx2xU'(x), . 1 ^/2X k U'(X k ) 



f dx 2xU'(x) , , 1 v-^ 



\2 

k 



I dx 2xU\x) , N 1 ^/ 2X k U'(X k ) -1 \ 

and insert them into the equations (3.33) and (3.43). This results in the following 
constraints, which are equivalent with the loop equations (3.33) and (3.47): 

< w « V-*) + «(*M-*)'> - ^ <^>+^>- 2 "<°» (3.49) 
/ da; 2xU'(x) , , . . 

3.4 The large iV Riemann surface 

In this section we shall take the large N limit of the loop equations for the orientifold 
with symmetric and antisymmetric matter, thus obtaining two polynomial constraints 
on the planar limit uq of the averaged resolvent 8 , which is the leading term in the large 
N expansion: 

(cu(z)) = u (z) + ^cu 1 + O(l/N 2 ) (3.50) 

These large N constraints on uj define a Riemann surface, similar to what happens in 
the much better studied case of one-matrix models. The main difference is that the 
Riemann surface governing our models is not hyperelliptic. Rather, it will turn out to 
be a triple cover of the z-plane, one of whose branches gives u> (z). 

8 These constraints also follow from the loop equations of the ^-quiver, which are given in Appendix 
A.2.1. 
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3.4.1 The loop equations at large iV 

The large N limit of the loop equations (3.33) and (3.47) gives: 

uj (z) 2 + uj (z)uj (-z) + M-z) 2 - V'(z)uj (z) - V'(-z)uj (-z) + f(z) + f(-z) = 
lo (z) 2 lo (-z) - g(z) + V\z){u {zf + f(z) -V'(z)u (z)) + (z^-z) = (3.51) 

It is convenient to introduce the shifted variable u(z) through the relations: 

cuq(z) = u(z) — t(z) , uo(—z) = u(—z) — t(—z) , (3.52) 

where 

t(z) = t( _ z) = -m-pnA (3 53) 

Under this translation, equations (3.51) become: 

u(z) 2 + u(z)u(—z) + u(— z) 2 = p(z) 

u(z) 2 u(-z)+u(z)u(-z) 2 = -q(z) , (3.54) 

where: 

p(z) = t{z) 2 + t(z)t(-z) + t{-z) 2 - f(z) - f(-z) 

q{z) = -t(z)t(-z) [t(z) + t{-z)\ + t(z)f(-z) + t(-z)f(z) - g{z) - g(-z) (3.55) 

Equations (3.54) are the Viete relations for the roots Uo(z) := u(z),Ui(z) := u(—z) and 
u 2 {z) := —uo(z) — Ui(z) of the cubic: 

2 

Y[(u(z) - Ui{z)) = u 3 - p(z)u - q{z) = . (3.56) 

i=0 

Therefore, the quantities u (z), u±(z) and u 2 (z) are the three branches of the affine 
curve (3.56), when the latter is viewed as a triple cover of the complex z-plane. 
Let us index the branches such that: 

ivo(z) = uq{z) — t(z) and uq(—z) = ui(z) — t(—z) . (3.57) 

It is clear from (3.55) that the functions p(z) and q(z) are even. Hence the curve (3.56) 
admits the automorphism: 

v(z,u) = (—z,u) , (3.58) 

which permutes the sheets u and u±, while stabilizing the third sheet. Using (3.55), 
our curve can be written: 

(u - t(z)) (u - t(-z)) (u + t(z) + t(-z)) 

+ [f( z ) + f(- z )] u + [ 2 g ev (z)-t(z)f(-z)-t(-z)f(z)} = ,(3.59) 
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where: 

9ev(z) := \ [g(z) + g(-z)\ = l - [g(z) + g(-z)\ (3.60) 
is the even part of g. For / = g = 0, the defining equation reduces to: 

(u-t (z))(u-t 1 (z))(u-t 2 (z))=0 , (3.61) 

where t (z) := t(z),ti(z) := t(—z),t 2 (z) = —t (z) — t\{z). Correspondingly, the 
branches are given by Uj = tj in this limit. The general curve (3.59) is a deforma- 
tion of (3.61), parameterized by the coefficients of / and g ev . 

The number of independent coefficients in / and g ev is constrained by the matrix 
model. Recall from the discussion of equations (3.37) and (3.45) that the polynomials 
f(z) and g(z) have degree d—1 if the polynomial part W(x) of the potential has degree 
d+1. Thus f(z) depends on d complex parameters while the even polynomial g ev 
depends on 5+1 parameters, where 

6= p*=i] . (3.62) 

In the large A" limit the function f(z) in (3.35) has the form f(z) = f(z) — (f>t-i/z, so 
it depends on the coefficients of f(z) as well as on 0, if the logarithmic term is present 
in the potential. Altogether, we have: 

(d + 5+1 ;t_i = a.) 
#( coefficients in /, g ev ) — < (3.63) 

[d + 5 + 2 ;t_i^0 b.) 

3.5 The Riemann surface in the absence of a logarithmic deformation 
3.5.1 General description and parameter count 

In the absence of a logarithmic deformation (t-i = 0), the large algebraic curve 
is given by (3.56) with V(z) = W(z). We shall start by counting its parameters 
and periods. For this, we must describe the branching structure of this curve and of 
its deformations. For a generic deformation (/, g), we shall find d + 5 + 1 independent 
periods, in agreement with the parameter count for g ev and / performed in the previous 
subsection. 

Let us start with the classical curve (3.61) and analyze its branching. Writing 

W(z) = Y / d m=0 tmZ m , we have: 

W' oM {z) := l -[W\z) - W'(-z)} =J2t2 k+ iZ 2k+1 = zv(z 2 ) , (3.64) 

fc=0 
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where we defined v(x) := J2i=o^+ix k . Consider the factorizations: 

W\z) = t d (z - 4) ...(z- 4) -> W'(-z) = t d {-l)\z - *D ...(*- *J) (3.65) 
with z~ := and 

v(x) = t 2 s+i{x - xi) . . . (x - xs) ■ (3.66) 
The second relation implies the factorization: 

W odd {z) = t 2 s+iz{z 2 - Xl ) . . . {z 2 - x s ) = t 2S+1 z(z - zt){z - zf) ...(z-zj){z- zj) , 

(3.67) 

where = rk^Jxj. 

Then the curve (3.61) has ordinary double points at those values of z where U(z) = 
tj(z) for some < % < j < 2. Let us assume that all roots zf,z~,z~^ and z~ are 
mutually distinct and nonzero (this is the generic case). Then: 

t (z) = t 2 (z) & W\z) = z e {z+ . . . z+} 

t^z) = t 2 (z) & W'(-z) = <£> z e {zi . . . z~} (3.68) 
t (z) = h(z) <£> W' odd {z) = & z e {^r . . . zj, 0, zt ... z+} . 

When turning on generic deformations f,g to reach the curve (3.56), all of these 
double points will split into cuts: 



z~l -> [a,i, bi] 
z+ -> [aj,bj] 
z -> [d 0: b ] 



h , ^ -> [-bi, -Oi] := I-i , i = l,...,d 
= Ij , 57 -> [-6j, -dj] := Lj , j = 1, . . . , 5 
= Jq , with a = — &o • (3.69) 



This splitting is symmetric since the allowed deformations must preserve the symmetry 
of our curve. We obtain 2d cuts of type ij and 25+1 cuts of type Ij. 

The deformed curve (3.56) has a multiple point at z = oo. It follows that an ap- 
propriate deformation of its normalization will have two supplementary cuts. Applying 
the Hurwitz formula for this normalized and deformed curve, we obtain its genus g: 

g = I(45 + 4d + 6) -3 + 1 = 2(d+5) + 1 . (3.70) 

This is one unit greater than the arithmetic genus of the curve (3.56). 

Choose cycles A i: Aj around the cuts ij, Ij as indicated in figure 1. Due to the Z 2 
invariant perturbation of the cuts and the induced action u*(X) = —A of (3.58) on the 
meromorphic differential 

\=J-udz, (3.71) 
2ixi 
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A-i 



A 



A 



A n 



An 



Figure 1: The cycles Ai and Aj. Solid lines lie on the sheet uq(z), while dashed lines lie on 
ui(z). Under the automorphism v the cycles transform as v{AA = A-i, v(Aj) 
u(A Q ) = -Aq. 



A_j and 



the periods integrals of this 'regularized' curve along A_i, A-j are minus the ones along 
Ai, Aj for % = 1 . . . d and j = 1 . . . 5, while J"^ A is invariant. Therefore, we can choose 
f A A and A (with i — 1 . . . d and j = ... 5) to be independent periods. Hence their 
number agrees with the parameter count for g ev and /. In (3.69) we assumed a generic 
situation, which can be justified from the double cover. In that case, the existence 
of 3d independent 3-cycles in the Calabi-Yau geometry has been established [36]. By 
the classical results of Tian and Todorov [43, 44], this is equivalent to 3d unobstructed 
complex structure parameters, which descend to the Riemann surface. The argument 
above can be viewed as establishing consistency of the Z2 projection on the parameters 
and periods of the Riemann surface. 

3.5.2 Physical meaning of the cuts 

We next discuss the interpretation of the cuts in terms of the planar eigenvalue distri- 
bution of the matrix model. Remember that uj (z) = J d\j^. From (3.57) we see that 
the cuts of u)q{z) coincide with those of Uq(z). The fact that Uq has branch cuts along 
Jj and Ij (with % = 1 . . . d and j = —5 . . .5) requires that p(A) be non- vanishing along 
each of these cuts. This implies that ujq(-z) will have cuts on the reflected loci J_j and 
I-j. To find the matrix model meaning of these cuts, consider the analytic function: 

k(z) := u (z) - u 2 (z) = 2u (z) + u (-z) = 2lu (z) + Uq{—z) - W'(z) = 

2 1 



dX'p(X') 



W'(z) 



z-X' z + A' 
where we used relations (3.52) and (3.53), which imply: 

2t(z)+t(-z) = -W\z) 

Similarly, we consider: 



(3.72) 



(3.73) 



t(z) := uq(z) 



ll {Z) - Ui(z) 

dX'p{X' 



u (z) - u {-z) = cu (z) - uJoi-z) - [W\z) - W\-z)\ 



+ 



X' z + X' 



- [W'(z) - W\-z)\ 



(3.74) 
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U 2 



A + ie 



A — ie 



h 



Uq 

Figure 2: The values of uq(X — iO) and U2^X + iO) agree along a cut Jj, since the latter 
connects the branches no and Ui. A similar argument holds for the other types of cuts. 

where again we used (3.52) and (3.53). 

The definition of the cuts /, and Ij implies (see figure 2): 



Uq(X ± zO) = «2(A T ^0) for A G /, with i = 1 . . . d 

Ui(X ± iO) = u 2 (X =F *0) for A G I_j with z = 1 . . . d 

«o(A ± ^0) = «i(A =F iO) for X E Ij with j = —S...S 

Therefore, we have: 

V(k(X)) = for A G Ii with i = l...d 

p( K (-A)) = for AG I-i with i = 1 . . . d 

V(t(X)) = for AG Ij with j = -6...8 . 

Writing these principal values using (3.72) and (3.74), we find: 

2 1 



dA'p(A') 
dX'p(X') 
dX'p{X') 



X 


- X' 


A 


+ A' 




2 




1 


X 


+ A' 


A 


-A' 




1 


+ — 


1 



A -A' A + A' 



W\X) , AG/,, 
-W'(-X) , AG/-,, 
W'(X) - W'(-X) , Xelj . 



(3.75) 



(3.76) 



(3.77) 
(3.78) 
(3.79) 



Notice that the second equation is equivalent with the first under the substitution 
A — > — X and therefore contains no new information. 
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Assuming that all cuts are disjoint from each other, the analytic function u 2 must 
be continuous across Ij (since the cuts of u 2 are Ij and /__,-). Therefore, we have the 
following relations for A e If 

u 2 (X + iQ) = u 2 (X - iQ) w (A + iQ) +u 1 (X + iQ) = u (X - iQ) +u l (X- iQ) ^==> 
u (X + iO) -uo(A-iO) = ui(A-iO)-ui(A + iO) . (3.80) 

Using tio(^) = 2) in the right hand side gives: 

p(A)= P (-A) , X e Ij , (3.81) 

where we also used the identity: 

u {X - iO) - u {X + iO) = 27rip{X) , XeR . (3.82) 

Thus p takes symmetric values along the cuts Ij and 7_j (in particular, p(A) is sym- 
metric along Jo)- 

Finally, we can identify the filling fractions of the matrix model: 

Si — / A , for % — 1 . . . d , 

Sj = 2 f x , for j = 1 ... 5 (3.83) 
S = [ 

J Ac 



X . 

'A 

The factors two in the second equation follows from the fact that Ij and support 
symmetric distributions of eigenvalues for j ^ 0. 

From the gauge theory point of view, we can interpret (3.77) as the quantum 
correction to the classical vacuum configuration (2.8) and (3.79) as the correction to 
(2.7). The different filling fractions correspond to the different types of vacua we found 
in Section 2.1. In particular, a fraction Sj with j ^ corresponds to a vacuum in 
which the final U(Nj) gauge group is embedded diagonally, while So corresponds to 
the vacuum with an orthogonal or symplectic unbroken gauge group. Of course, up to 
now we imposed the strict large iV limit and therefore the curve can not distinguish 
between the two types of orientifolds. In section 5 we will argue that this difference 
can be accounted for by including a logarithmic term in the potential V of the model. 
Let us therefore study this situation next. 
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3.6 Riemann surface in the presence of a logarithmic interaction 

In this subsection, we analyze the Riemann surface (3.56) in the presence of a logarith- 
mic term t-ilnz in the potential V. We shall show that turning on such interactions 
produces supplementary cuts, which correspond to eigenvalues accumulating along new 
loci. The purpose of the present section is to identify these novel cuts for small values 
of t-i — this will important in our discussion of the first order loop equations in section 
5. 

To understand the generic situation, first note that p(z) and q(z) acquire poles 
at the origin for t-i^O. To display this, we decompose t(z) = t(z) — — , with 
t(z) := ~ 2W ( z ^ w a polynomial. Further, we define polynomials p(z) and q(z) 
by substituting /, g, t in (3.55). This results in: 

M=P(z)~(i(z)-i(-z)) + t 4t 

A, jC 

q(z) = q(z) + t -^{f(z)-f(-z) + [i{z)-t{-z)\[<t>-t{z)-i{-z)\) (3.84) 

z 

where we used (3.36) and (3.45) and absorbed the s dependence in t-i. To write (3.56) 
in polynomial form, we multiply the equation by z 3 and define a new variable y := zu: 

y Z _ p( z )y _ Z Q( Z ) = ( y _ y ( Z ))(y _ y^y _ y^ z)) = Q , ( 3 .85) 

where yi(z) = -y (-z), y 2 {z) = -y (z) - yi(z) and P(z) = z 2 p(z) and Q(z) := z 2 q(z) 
are even polynomials. This curve admits the symmetry: 

( Z ,y)^(- Z ,-y) . (3.86) 

We shall view (3.85) as a deformation of the curve: 

f - z 2 p(z)y - z 3 q(z) = (y - y (z))(y - yi (z))(y - y 2 (z)) = , (3.87) 
which is obtained from: 

u 3 — p(z)u — q(z) — (u — u (z))(u — ui(z))(u — u 2 (z)) = , (3.88) 

where u 2 (z) = —u (z) — Ui(z) and Ui(z) = u (—z), by performing the birational trans- 
formation y = zu. 

We are interested in the generic behavior of (3.85) as t_i approaches zero. In this 
limit, the curve degenerates to (3.87), whose first two branches yo and y~\ are tangent 
in a cusp at y = z = 0, where they also meet the third branch y 2 (figure 3). 
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Figure 3: Branching of the curve (3.85) at the origin for t-i = 0. 

The geometry of (3.85) for a small value of t-i is described in figure 4. The loga- 
rithmic deformation of the model generates four new branch points of the birationally 
transformed curve, which are indexed by z± and z±. These are symmetric with respect 
to the origin, namely z_ = — z + and 2- = —z + . 











\Z-/ I 
















J / vo ^ 









Figure 4: Deformed curve after addition of a logarithmic term in the matrix model potential, 
in the birationally equivalent coordinates (y,z). 

From this geometry, one finds that the cut Jo splits into two new cuts Iq = [z + , bo], 
and Iq = [So, z~], which are distributed symmetrically with respect to the origin. Both 
of these cuts connect the branches and 1. One also finds two new cuts I = (— oo, z + ) 
and —Jo = (-2-, +oo), which connect the pairs of branches (0, 2) and (1, 2) respectively. 
As in [5], we shall 'regularize' our curve by choosing a large A > and replacing these 
cuts with I = (—A, z + ) and —I = (z_, +A). 
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These cuts produce new A-type cycles A , Aq (and associated B-cycles), which 
project onto the z-plane to the curves 7o,7o = , see figure 5. The homology classes Aq 
are interchanged by the symmetry z —>■ —z, which also interchanges the homology class 
of Aq with that of — Aq. 

The cuts I and Iq give the A-type periods: 



' 70 



To 



dz 






Ja 2™ ^ ' 


dz . . 


/ dz 


2^ oW = 





(3.89) 



The cut —Iq gives the opposite of the Ao-period. The Z2 symmetry also implies: 



To 



dz 
2tt7 



uj (z) 



dz 
2txI 



Uq(z) 



(3.90) 
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an I z. 
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• x 



h 
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7 



z+ 1+ b a 



z = 



Figure 5: The new cuts of uq generated by a logarithmic term in the potential (viewed as 
cuts of the birationally transformed curve (3.85)). The cut —Iq corresponds to branching of 
Ui with U2, and consequently it is not shown. The figure also indicates the new cycles Aq, 
Aq on the Riemann surface. Solid cycles lie on the sheet uq(z), while the dashed cycle lies 
on u\(z). The curves 70, 7^ mentioned in the text are the projections of Aq and onto the 
z-plane. 

The presence of new cuts for the log-deformed model shows that eigenvalues can 
accumulate along I and Iq in the large N limit (no eigenvalues accumulate on — J since 
this is not a cut of Uq) . Accordingly, the microcanonical ensemble of the log-deformed 
model involves the supplementary constraint: 



10 



dz 
2ni 



(u(z)) = Sq 



(3.91) 
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while the constraint §~ -^(uj(z)) = S of the undeformed model is replaced by: 



This becomes part of conditions (3.27) for the case t-i 7^ 0. For a consistent limit 
as t-i — > 0, one takes So to be a function of i_i which tends to zero in that limit. Note 
also that with the additional constraint (3.91) we have to fix d + 5 + 2 periods in the 
log-deformed model, in agreement with the parameter count in (3.63). 

4. Chemical potentials at large N 

In this section, we study the planar limit of our models for the case i_i = 0. Our 
main purpose is to give a proof 9 of an appropriate set of special geometry relations. 
The path we shall follow is based on our construction of the micro canonical ensemble, 
combined with a modification of an argument due to [45]. Namely, we shall show 
that the chemical potentials \i reduce in the large N limit to certain B-type periods 
II of our Riemann surface. Then special geometry conditions of the type found in [1] 
follow from the standard equation (3.26) which expresses the chemical potentials as 
derivatives of the microcanonical partition function with respect to the filling fractions. 
In particular, equation (3.26) is the appropriate finite N generalization of the special 
geometry constraints. We also show that Whitham-type relations of the type found in 
[45] arise as the planar limit of certain finite N equations which follow naturally in the 
microcanonical ensemble. For simplicity, we assume t-i = for most of the present 
section. The generalization to the case t-i 7^ is entirely obvious, but notationally 
tedious. 

Let us hence assume t-i = 0. Since p(A) in our models develops a nonzero value 
along the cuts ij and Ij with % = 1 . . . d and j = —5 . . .5, we have d + 25 + 1 filling 
fractions. To construct the microcanonical ensemble, we pick intervals A$ (i — 1 . . . d) 
and Aj (j = —5 ... 5), such that A_j = — Aj for all j. In the large N limit, we shall 
assume that the intervals Aj and Aj contain the cuts ij and Ij. The symmetry property 

9 The original paper [1] gives a beautiful intuitive argument for the existence of such a relation, 
without providing a rigorous proof (what is missing is to show that the Riemann surface B-periods 
indeed equal appropriately defined chemical potentials). A derivation of this relation was later given 
in [45, 46, 47] (using older results of [48]), though a clear construction of the microcanonical ensemble 
(which is implicit in that argument) was not given there. In the present paper, we are dealing with 
the more complicated case of a non-hyperelliptic Riemann surface, which underscores the usefulness 
of having a clear proof of such relations. 




(3.92) 
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(3.81) shows that the filling fractions of Ij and 7_j are equal, and we shall take So . . . S$ 
to be the independent quantities: 

/ d\p (X) = Si for % — \...d 

Jii 

J_d\p (X) = ^S j for j — 1...5 (4.1) 
/ dXpo(X) = S ■ 

Then the d + 8 + 1 constraints (4.1) completely fix the d + 8 + 1 deformations 
of (3.56) encoded by the polynomials / and g ev . In the (grand) canonical ensemble, 
this is implemented by introducing chemical potentials pi and jlj with i — 1 . . . d and 
j=0...5. 

The planar (grand) canonical generating function reads: 

T (t,p) = j dXW(X) P o(X)- 1 - J dxj dX'K,{\X')po{X)po{X') 

d 8 

+ £^3 + 5^/2^ (4.2) 

i=l i=o 

where 

/C(A,A') = 2 In |A — A'| -ln|A + A'| . (4.3) 
This gives the planar limit of the microcanonical generating function: 

F (t,S) = \Jd^J dX'IC(X,X')p (X) P o(X') - J dXW(X) P o(X) , (4.4) 

with the constraints (4.1). 

Observation Remembering relation (3.81), one finds the functional derivatives: 



8p(X) 
8F 



= J dX')C(X, X')po(X') - W(X) , Xel:= U t=1 ... d k 
= J dX'iC{X,X')po{X')-W{X)-W{-X) , \eI:=U j= -8..Jj , (4.5) 



8p(X) 
where: 

i<Z(X, X') = /C(A, A') + /C(-A, A') = In | A + A'| + In |A - A'| . (4.6) 

Note that equating (4.5) to constants and further differentiating with respect to A gives 
the planar equations of motion (3.77,3.78,3.79). This observation can be viewed as an 
intuitive justification for the rigorous procedure discussed below. 
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4.1 The primitive of k along the real axis 

To extract the large N chemical potentials, we shall be interested in the 'restriction' of 
the function k of subsection 3.5.2 along the real axis, which we define by: 

Kp(A) = i [k(A + iO) + k(A - iO)} (4.7) 

for any real A. If A is a real value lying outside the union of /«, then k p (X) equals k(X), 
the quantity obtained by substituting A for z in (3.72). Taking the principal value of 
(3.72) along the real axis gives: 



Kp(A):= J dX' Po {X')K(X,X')-W'(X) , (4.8) 
where K{z, z') := — Consider now the function : R — > C defined through: 

0(A) := J dX']C(X,X')p (X')-W(X) . (4.9) 
Noticing that X (A, A') = ^/C(A, A') shows that is a primitive of /Cpi 

As discussed in the previous section, the equations of motion (3.28) amount to the 
requirement that k p vanishes along each of the intervals If. 

k p (X) = for \ el = \jf =1 l i . (4.11) 

This means that is constant along each of these intervals: 

0(A) = & = constant for X E h . (4-12) 

Then the jump in the value of k p between consecutive cuts can be obtained by inte- 
grating (4.10) between their endpoints: 

6+1-6= dXK(X) . (4.13) 

Jb, 

4.2 The primitive of r along the real axis 

An entirely similar discussion can be given for the function t(z), whose 'restriction' to 
the real axis is: 

T P (A) = X - [r(A + tO) + r(A - iO)} = J dX' Po (X')K(X, A') - [W'(z) - W'(z')} , (4.14) 
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where K(z, z') : 



Z — Z' 



1 



7 Its primitive 



1 



d\')C(X, A')po(A') - W(X) - W(X') 



(4.15) 



can be integrated in between the cuts Ij to give: 




dXr(X) , 



(4.16) 



where ^ is the constant value of ip along Ij (as required by the second planar equation 
of motion): 



4.3 The large N chemical potentials 

Differentiating (4.4) with respect to Si for some i < d and using relation (4.9) gives: 



where we took the dependent filling fraction to be Sd- 

To arrive at the last equality, we used equation (4.12) and the constraints (4.1) and 
(3.25). Relation (4.18) shows that the chemical potentials associated with Ij coincide 
with the quantities £j in the planar limit, up to a common additive constant. Using 
relation (4.13), we obtain: 



A similar argument for the differential with respect to Sj (using relation (4.15) and 
the last two constraints in (4.1)) gives jif^ = |j and: 



r p (X) = for X E I = \J j= -s...sIj 



(4.17) 




(4.18) 




(4.19) 




(4.20) 



4.4 Geometric expression for the large iV chemical potentials 



Together with (3.72), equation (4.19) gives: 




(4.21) 
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where Bi are cycles on the large N Riemann surface chosen as explained in figure 6. 
Similarly, we find: 

fif - = C +1 dz[u x {z) - u (z)} = I dzu(z) , (4.22) 
Jlj Jb 3 

with Bj chosen as in figure 6. 



'd-2 



A 



d-2 



B, 



d-3 



h-1 



A 



d-l 



A d 



B, 



d-2 



B 



d-l 



B d (A) 



A 



Bs-i 



B S (A) 



Figure 6: Choice of B-cycles on the large N Riemann surface. We only indicate a few 
cuts close to the point A. The cycles Bi,Bj and Ba(A), Bg(A) are defined such that, when 
crossing a cut going upwards along these cycles, one moves from the branch uq to the branch 
u\ (drawn dashed) or U2 (dotted) respectively. With the orientation of the Riemann surface 
induced by its complex structure, and with the indexing a of cuts explained the text, this 
implies the intersections A a r\B a = +1 for the re- indexed cycles (note that the cycles A a and 
B a intersect in a single point, which lies on the branch uq). 

Consider the cycles Bi = Ylt=l fo r all z = 1 . . . d — 1. Then (4.21) implies: 

lif = {if + / dzu{z) for i = 1 . . . d - 1 . (4.23) 

JBi 

The quantity ^ is undetermined and can be fixed arbitrarily. Following [1], we take 
rd° = §B d dzu(z), where A is a point close to +oo and Bd(A) is a path shown in figure 
6. Defining Bi(A) = Bi + B^(A) for all i — 1 . . . d — 1, equation (4.23) gives: 

= Tli for i = l...d , (4.24) 

with: 

n, := / dzu(z) . (4.25) 

JBi 

A similar construction for the cuts Ij gives: 

=tl j for j = ... 5 , (4.26) 
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where: 

TLj := / dzu(z) , (4.27) 
with cycles Bj chosen analogously (figure 6): 

5-1 

Bj = B 5 (A) + Y,Bk ■ (4.28) 

k=j 

Relations (4.24) and (4.26) show that the chemical potentials \x are the finite A" ana- 
logues of the periods II. 

Also remember that the filling fractions can be expressed as periods of udz over 
the cycles A a of figure 1: 

S. = j |W„ W = j ±u o(z ) = j ^uiz) . (4.29) 

In the second equality, we used relation (3.52) and the fact that t(z) is a polynomial. 

Defining q :— d + 8, let us index the cuts by ct = — q ... g in the order in which the 
intervals I iy Ij appear from left to right. We also let A a and B a respectively B a be the 
corresponding cycles A iy Aj and B iy Bj respectively B^ Bj indexed in this order. It is 
then clear from figure 6 that A a DB a = — A a C\B a -i = +1 and A a C\Bp = if (3 ^ a, a— I. 
This gives A a fl Bp = —Bp fl A a — S a p. Since we also have A a H Ap — B a n Bp — 0, it 
follows that A a , Bp have canonical intersection form. Thus we find a canonical system 
of cycles A a , B a with a — — q . . . q : 

A a (~)Bp = -B a nAp = 5 a> p , A a nAp = B a nBp = for all a = -g . . . q . (4.30) 

4.5 Whitham-type relations 

Let us now show how a finite A" version of a set of Whitham-type constraints (similar 
to those found in [45] for the case of one-matrix models) can be extracted from the 
microcanonical ensemble. For the rest of this subsection, we shall allow for a logarithmic 
term in the potential — as we shall see in a moment, this has interesting effects on our 
Whitham-type relations. 

Differentiating (3.24) with respect to t rn gives: 

-F(t,S) = -— Ht,^ a =, a (t,s) = lm , (4.31) 



dt m dt r 



where we defined the averaged loop operators: l rn := —-^(trd tm V(M)) = — at(t^+i) (tr(M m+1 )} 
for m > and l_i := — j^{tr{logM)) for m = —1. Equations (4.31) supplement rela- 
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tions (3.26). We have: 



_ z = f^Jd\m)^ = £- 1 f J £z^(u,(z)) ifm>0 
\jdX(p(X))\ogX^^^(u;(z))\nz if m = -1 

where 7 is a counterclockwise contour that encircles all eigenvalues. Combined with 
(4.32), equations (4.31) give a finite N version of a set of Whitham-like constraints; 
these reduce to standard Whitham conditions in the planar limit and for the case 
f_i = 0. 

4.5.1 Geometric form of the large N Whitham-type constraints 

For m > 0, one can deform the contour 7 toward infinity, thus picking up contributions 
from z = 00. This gives: 



lm = —ReS(; = o 

m + 1 



HVO) 

c 



m+3 



for m > , (4.33) 



where ( = 1/z. 

In the large N limit, one has (w(z)) = u>o(z), with ouo(z) a branch of the Riemann 
surface (3.56). In this case, equations (4.33) and the last equation in (4.32) determine 
the derivatives |p in terms of the coefficients of V and /, g; these are Whitham-type 
conditions analogous to those of [45]. As shown above, the derivatives are also 
determined by the Riemann surface. Together with the Whitham relations, this allows 
one to determine the planar generating function Fo(t, S) up to an additive constant 10 , 
given the planar limits l^(t, S) and U(t, S). 



5. First order analysis of the loop equations 

Consider the first order large N expansion of the microcanonical generating function 
F(t, S) of the log-deformed model: 

F(t,S) = F (t,S) + jfF&S) + 0(N~ 2 ) . (5.1) 

The aim of this section is to show that the contribution Fi can be obtained by differ- 
entiation with respect to the parameter t-\ in the potential (3.11): 

= S z dJ £f- + »(*-!. *) ■ (5-2) 

10 We stress that equations (3.26) only determine Fo(t, S) up to the addition of an arbitrary function 
of t m . One needs the Whitham constraints in order to fix the t-depcndcnce of F . 
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5.1 The resolvent to order 1/N 

To establish this, we first note that the multi-point correlators of the resolvent uj differ 
from the corresponding product of one point functions only at order 1/N 2 : 

{u{z) 2 ) = (u J {z)) 2 + 0{l/N 2 ) , (,(,) 2 ,(-.)) = <^)>V(-^)>+0(l/iV 2 ) (5.3) 

and so forth. Consider the 1/iV-expansion: 

{u>(z))=u (z) + ±u> 1 (z) + O(l/N*) . (5.4) 

Inserting this and (5.3) into the integral formulation (3.49) of the loop equations we 
obtain the planar relations: 

dx 2xV (^) / \ /\2 / \ / \ / \2 

wo{x) = UJ {Z) + uj (z)ujq(-z) + LJq(-Z) , 



2m z 2 — x 2 
dx 2xV'(x 



fiy 2ni z 2 — x 2 



u> (x)u> (-x) = u (z) u (-z) + uj (-z) u (z) (5.5) 



and the 0(1/N) constraints: 

dx 2xV'(x) I dx u (x) 



f dx 2xV'(x) . , / 



7 2m z 2 — x 2 



= [2u (z) + u (-z)]ui(z) + [2u Q (-z) + u (z)]ui(-z) , (5.6) 
dx2xV'(x) T f dx u (x)uq(—x) 



2m z 2 — x 2 ./I 2ixi z 2 — x 2 

= [2ujq(z)ujq(-z) + u)q(-z) 2 ]uji(z) + [2u (-z)u (z) + lv (z) 2 ]ui(-z) . 

On the other hand, the filling fraction conditions give: 
and 

£"* = ° (i = '"' i) • |^ = {3 = 1 - S) • i^ = ' (5 ' 8) 

These include the constraints for the log-deformed model arising for the cuts shown in 
figure 5. 

Suppose that cu (z) is a solution to (5.5) and (5.7) for the potential (3.11). By 
differentiating the planar loop equations (5.5) as well as (5.7) with respect to t-\ one 
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can check that lo\ := fj^- is a particular solution of both (5.6) and (5.8). To show that 
this is the unique solution, it suffices to check that the homogeneous system obtained 
from (5.6) by dropping the terms s ^ |^£§r and s ^ —^z^ 1 from the two left 
hand sides admits only the trivial solution when supplemented by the constraints (5.8). 
To establish this, notice that the homogeneous system associated with (5.6) results from 
the planar equations (5.5) if one performs infinitesimal variations of u\ = 5ujq, such 
that the variations are independent of the coefficients £_i . . .t^ of the potential (3.11). 
Such variations of ujq arise by changing the filling fraction parameters S a , which give a 
parameterization of that part of the moduli space of our Riemann surface which results 
by varying the quantities /, g of Subsection 3.3. Therefore, the solution space to the 
homogeneous system associated with (5.6) is the tangent space T to this part of the 
moduli space of our surface. Since the corresponding moduli space is parameterized 
by S a , the dimension of this solution space equals d + S + 2. On the other hand, 
relations (5.8) give d + 5 + 2 linearly independent equations for such solutions. These 
equations select the vanishing solution of the homogeneous system associated with (5.6). 
Therefore, the homogeneous system obtained from (5.6) has no nontrivial solutions 
when supplemented by the constraints (5.8). This shows that the solution of (5.6) and 
(5.8) is uniquely determined and equals 

We conclude that the solution of the loop equations has the form: 

(lo(z)} = Uo (z) + + 0(1/N 2 ) , (5.9) 

where: 

sdu (z) 

= • (5.10) 

5.2 The 1/N correction to the microcanonical generating function 

Next we show that a similar relation holds for the 1/iV-expansion of F. To simplify 
notation we label all filling fractions S, S by S, keeping the convention that So denotes 
the filling fraction of Iq in figure 5. Equation (5.10) was derived under the assumption 
that the nonsingular part W of our potential contains a finite number of terms d. Since 
(5.10) must hold for any value of d, it is clear that it also holds if we formally allow 
W to contain an infinity of terms. Therefore, let us now take W(z) = J2m=o ~^+i zm+1 ■ 
Defining the loop insertion operator dV ^^ = Xlm=o 1^^W~ ( see e -S- [49])> we then 
have (m + 1) Jg- = ±(trM m+1 ), for all m > 0, which gives the standard relation: 
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upon expanding u(z) = ^(tr^) = ^ Em=o P^t^PT^- Here F is the micro- 
canonical generating function. Expanding F = X^>o jjjFji this gives tuj = — dV ^^ Fj . 
Therefore, relation (5.10) implies: 

Fi(t-i, ...t oo ,S ...S oo ) = s _ dF o(t-i---too,S ...S 00 ) + ^ So ^ Soo) (512) 

where \1/ is a universal function of t-i, So . . . which does not depend on (tj)j>o- Since 
this function is completely independent of to . . .too, it knows nothing about which of 
these coefficients are zero. To determine it, let us consider a potential (3.11) with 
W = 0; this amounts to setting tj = for all j > 0. In this case, the associated 
Riemann surface (3.85) has a single independent period, namely So- It follows that the 
function \1> does not depend on the variables Si... S^. Thus one has \I> = Sq). 
Returning to equation (5.12), we now take W to be of the form (3.11) and obtain: 

Fl (t^, t ...t d ,S ... S d+S+1 ) = S - dFo{t ^ to • • • • • Sd+5+l) + *(t- u So) . (5.13) 

Since we are ultimately interested in the case t-i = 0, let us set S'o = S (t-i) with 
So(0) = in equation (5.13). Then the limit t-i — > gives: 



iq(0, t • • • td, 0, Si . . . S d+ s+i) = - 



dF (t-i,tp . . . td, S . . . S d+ s+i) 



+^(0,0) . 

t_i=0,S =0 

(5.14) 

Therefore, the microcanonical generating function of the model with t_ x = satisfies: 

s dF 



F l 



2dt-i 



t_i=0 



+ constant . (5.15) 



This relation is reminiscent (though somewhat different in character) of similar equa- 
tions satisfied by the (microcanonical) generating function of SO(N) and Sp(N/2) 
one-matrix models [22, 28]. 



6. Comparison with the field theory description 

In this section we want to compute the leading terms in the effective superpotential 
from our matrix model. The strategy to do so will be the following. We will take the 
superpotential of the gauge theory in which all the fields are complex. So to start with, 
there is no hermiticity constraint on M and also Q and Q are unrelated. The gauge 
symmetry will be given by the complexified gauge group just as in the discussion of 
the classical moduli space of the gauge theory. After imposing an appropriate gauge 
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fixing a la BRST we then choose a suitable real section in the configuration space and 
perform the path integral along this real section. This is completely analogous to the 
situation in [7, 18] where in a two-cut vacuum of the one-matrix model one was forced 
to expand around Hermitian and ant i- Hermit ian matrices respectively. This somewhat 
ad hoc procedure has been justified recently in [20]. 

According to the DV-conjecture the leading logarithmic terms in the superpoten- 
tial can be computed from the "non-perturbative" part of the matrix-model partition 
function. Actually the notion of non-perturbative should be understood here in the 
sense that we compute the one-loop approximation to the partition function, taking all 
couplings of terms involving the product of more than two matrices to zero. In other 
words, we will expand the matrix model around the background only up to terms that 
are bilinear in the fluctuating fields. In order to keep things simple but also capture 
the essential physics we choose to expand around the vacuum with vevs given by 

(M) = diag(0jv , al Nl , bl N2 , -bl N2 ) ; (6.1) 



and 



(Q) 



V 







1 N2 
s!n 2 



(Q) 



f-sW'(0)E s 



\ 



\ 









-W'(b)l N2 



sW'(b)l N2 




(6.2) 

The N x N matrix E s is defined as in Section 2. The constants a and b are solutions 
to the equations W'(a) = and W'(b) = W'(—b). This vacuum is the simplest generic 
one in the sense that each of the different types of vacua we found in the analysis of 
the classical moduli space appears precisely once. 



6.1 Fixing the gauge a la BRST 

The gauge transformation are 

5M = i[A,M], 5Q = i(AQ + QA T ) 5Q = -i(A T Q + QA) (6.3) 

Since we are interested only in the quadratic part of the action we can also linearize 
the gauge transformations around the background. Furthermore from now on we will 
decompose all matrices in blocks of sizes N x Nj with i,j e 0---3 and denote the 
fluctuations of the matrices around the vevs with lower case letter, i.e. M = (M) + to, 
Q — (Q) + 1 an d Q — (Q) + Q- The linearized gauge transformations then take the 
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form 



Srriij = - (M) ii )A ij 

Sqij = i(A ik (Q)kj + (Q)ik(A jk ) T ) 

hij = -i{(^ki) T (Q)kj + (Q)ifcAjy) 

From these transformation we see that we can choose the gauge 

rriij = ; for j ^ i , 
Q23 = goo = . 



(6.4) 
(6.5) 
(6.6) 



(6.7) 



The gauge fixing is now implemented in a standard fashion by promoting the gauge 
transformations to BRST transformations, i.e. replacing the gauge parameters by ghosts. 
We also need to introduce antighosts C and Lagrange multiplier fields B. The gauge 
fixing Lagrangian is given by 



S gf = sTr (y^^Cijiriij) + C q q 23 + C 00 <7oo^ 



(6.8) 



The BRST transformations for the antighosts and Lagrange multipliers are 



sC q = B q , sB q 

s Coo = -^00 ! s -^00 

Explicitly the gauge fixing action is 

S 9f 



= 0, 

o, 

= 0. 



(6.9) 



Tr [ Bijmij + i((M)u - W 3j )C l3 C l3 + B q q 23 - tC q (C 22 + Cj 3 ) + 

i¥=3 



Booqoo — iCoo(CooE s + E s Coo) 
The the residual gauge group is 

G s = l/(JVi) x U(N 2 ) x 



SO(N ) s = +1 
Sp(f) s = -1 ■ 



(6.10) 



(6.11) 



Notice that the U(N 2 ) group is diagonally embedded in U(N 2 ) x U(N 2 ) such that the 
rank of the original U(N) gauge group is N = Ni + 2N 2 + N . According to this 
structure of the gauge group only the combination C 22 + C 33 appears in the gauge fixed 
action. Similarly only the ghosts components of C o that do not correspond to SO/Sp 
residual gauge transformations propagate. 
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6.2 Gaussian approximation of partition function and microcanonical gen- 
erating function 

Expanding the action around the background up to terms that are bilinear in the fields 
gives 



S (2) = N ( AW(a) + N 2 (W(b) + W(-b)) + N o W(0) + W"(( M >-) Tr (™«) 

\ i=0 

3 

]((M) U + (M)jj)Tr (<!,,<!„) + ^(M)^Tr + 



+ 



E 

i>j i=0 



+Tr (q 32 (m 22 + m^)) + -Tr (g 00 (m 00 E s + E s m^ )) . (6.12) 

In this expression we used already that q 2 3 = 0, goo — and m^- = for i ^ j. Notice 
that this gauge fixing is also important for another reason. Even without our gauge 
choice the blocks g 2 3, loo do not contribute to S^ 2 \ This is because {M) 22 = b = —(M) 22 
and (M)oo = 0. Fixing the gauge in the way we did eliminates this problem. Moreover, 
the fields q 23 and goo ac t now as Lagrange multipliers implementing the constraints 

m 22 = —mj 3 and m oE s + E s mQ = . (6.13) 

This restricts the fluctuations of M precisely in such a way as to account for a field in 
the adjoint representation of the diagonally embedded U(N 2 ) gauge group and another 
field in the adjoint representation of SO(N ) or Sp(^) respectively. 

In the definition of the matrix model in section 2 we had to implement also a 
reality constraint of the form M — = 2iel and Q = —iQ}. However, for the 
present calculation it is more convenient to choose different sections to ensure that all 
exponentials are decaying. For instance if a sits at a maximum of W(x) we will choose 
mn ant i- Hermit ian. Similarly, we choose q\j = qjj or q^ = —qjj depending on the sign of 
the coefficient of q^qij in S^ 2 ). For the ghosts we choose either Re (CV,) = Re (CV,) = 
or Im (Cij) = Im (Qj) = 0. 

In addition we can not relate all blocks of q to blocks of q by a conjugation. De- 
manding for example q 32 = q 2S would automatically set q 32 to zero because of the gauge 
fixing. We do need these fields however in order to implement the constraints (6.13) 
on the fluctuations of M. Therefore we will demand q 32 imaginary and goo imaginary 
symmetric or antisymmetric depending on s. After integrating out the Lagrange mul- 
tipliers Bij, B q , B 0Q , q 32 and goo the partition function of the gauge fixed matrix model 
is given by 



vo\(G s 



7q~T j W_ dm a II d( to dC d0 e (5(2)+Sfl/) (6-14) 
^ s) J ie{o,i,2} iy . I (23) 

(00) 



*{ 
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Performing the Gaussian integrations gives 



where 



J M 



^ — TTTTT^M Z ghosts % 'matter , (6.15) 

vol(G s ) 

/ 2tt \ 4 f 2tt \ 4 f 2tt \ ^(No-s) 

\W"{a)J \W"{b) + W"{-b)J \W"{0)) 
ghosts = (a 2 -b 2 f N ^a 2N ^{2b) 2N h m ^ 

Ka^) U Kb) ■ (6 - 16) 

The microcanonical generating function is defined through N 2 F = — log(Z). In order 
to compute it we need the logarithms of the volume of the residual gauge group. These 
are given by 

N 2 f 3 \ 

log(vol(C/(A0)) = — " g " log ( 27r ) ) + ( iY °) > 



Z, 



log(vol(/J s (iV))) = ~ (\og{N) - 1 - log(27r)) + 



sj (log(iV) - 1 - log(4) - log(27r)) + 0(N°) , (6.17) 



where H s is SO(N) for s = 1 and Sp(^) for s — — 1. According the the DV conjecture 
we introduce now the filling fractions Si — ^ and define the masses m(a) = W"(a), 
m(b) = (W"(b)+W"(-b)) and m(0) = W"(0). The 1/N expansion of the microcanon- 
ical generating function we denote by 

F = N 2 (f° + jyF 1 + 0(1/N 2 )^j . (6.18) 

The complete microcanonical generating function is 

F = F M + F ghosts + F matter (6.19) 



where 



Fm — 



Y log {^r—J -T log ~T log l^T7J + 

+ iV 4 8 Vem(0)V' 

2 

Matter = "^2 log (^^) " (^0 + Sf/2 - ^S 1 ) log Q - (6.20) 

-(2S 2 S + S 2 2 + sS 2 )\og(^y 
F gh0 sts = ~2S 1 S 2 log(a 2 - b 2 ) - 2S 1 S log(a) - 2S 2 2 log(26) - 4^0 log(6) . 
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6.3 The superpotential in the Gaussian approximation 



The effective superpotential of the gauge theory can be computed from these expressions 

as 



The subleading 1/N term in the microcanonical generating function stems form graphs 
with topology MP 2 . These are present in our theory because of the symmetric/ antisymmetric 
Q matrices. It has been argued in [21, 22, 28] that the MP 2 contribution to F enters the 
superpotential with a an additional factor 4. This argument has originally been derived 
in the context of theories with orthogonal and symplectic gauge groups. However, it re- 
ally only depends on the topology of the Feynman graphs and therefore should go over 
to our case without essential changes. In any case we will see shortly that (6.21) gives 
the correct Veneziano-Yankielowicz form of the superpotential. The final ingredient is 
a. It is given by the coupling of the original U (N) gauge theory 



The factor iV — 2s is the one-loop beta function coefficient of the original gauge theory. 
The scale /i is the scale where the potential W(4>) becomes relevant. For notational 
simplicity we will set this scale to one. Alternatively we could have introduced a 
scale in the matrix integral such that the Si obtain the correct dimension three as is 
appropriate for the interpretation as gaugino condensate. The ctj are then given by 
«o,i = ol and a 2 = 2a, where the last factor two stems from the diagonal embedding of 
the gauge group U(N 2 ). Collecting now all the terms we find the effective low energy 
superpotential 




(6.21) 




(6.22) 



W eff = Si log 



+5i log 




A' 



3JVi 
low,l 



+ S 2 log 



+ S 2 log 





\ So, / ge^ \ _ 



■N+2s 



7T 



2(N-2s)-2N 2 



The relation between the low- and high energy scales is 




AN 2 A 2 ( Ar - 2s ) 



(6.24) 



43 



Thus the effective superpotential coincides with the one obtained in section 2 up to the 
linear terms in the SV Of course these linear terms can not be determined by naive 
threshold matching. 
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A. Comparison with the 0(n) and matrix models 
A.l Comparison with the 0(n) model 

Let us compare our models with symmetric and antisymmetric matter with the 0(n) 
model analyzed in [42] . The difference is the (M, Q) interaction term, which reads in 
the 0{n) model 



where Q\ = Qi are Hermitian. After diagonalizing M in the relevant case for compar- 
ison 11 n = 1, the Gaussian integration over the real Qu leads in our conventions to 
s = 0. To summarize one has s = 1 for Qu complex, s = for Qu real and s = — 1 
for Qu = and this affects only the definition of U(x) (3.12) while the rest of the 
derivation of the loop equations remains the same. The quadratic loop equation (3.33) 
agrees in the large A-limit with the one in [42]. A cubic loop equation does not appear 
in [42] explicitly. The second order loop equations generalize very directly to general 
n. The third order loop equation should become an n + 1 order loop equation. From 
the gauge theory point of view the n = 1 model corresponds to adding matter in the 
adjoint representation, a spectrum which emerges e.g. from an Af = 4 theory. Adding 
more matter n > 1 leads to non asymptotically free cases. 

A. 2 Comparison with the A 2 quiver theory 

Let us next compare the situation to the ADE quiver theories, in particular to A 2 , 
which is the double cover of the models discussed so far. The ADE matrix models where 

11 For a special cubic V(x) this corresponds to Ising model on a random lattice [50]. 



n 




(A.l) 



i=i 
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discussed in [38, 39, 41] and further considered in the context of the DV conjecture in 
[2, 51, 52]. The action of the ADE quiver matrix models is: 



Sade = ^trV%(M«) + SijtiQ {ij) M^Q Ui) . (A.2) 

i=l i,j=l 

Here one has an x matrix for every node in the Dynkin diagram of the 
ADE Lie algebra and the Q(ij) are x matrices transforming in the (N^\ N^) 
representation and fulfilling = Q(j%)- For a given ordering of the nodes Sy = 

— Sji = 1 if node i and node j are linked in the diagram and s^- = otherwise. After 
diagonalizing all the Gaussians in the Q-fields appear with different signs. Let us 
assume = 1 and Qji = iQ\y Then we have to shift the eigenvalues into the 
upper half-plane whereas the M^) eigenvalues have to be shifted into the lower half- 
plane ! For Sij = — 1 we can still shift the eigenvalues into the upper half-plane 
but we then have to demand Qji = —iQ\y After choosing these real slices through the 
matrix configuration space the Gaussian integrations over all Q-fields can be performed. 

V T T 

N 2 S ADE = JV^E ^O^VEEMA^-A^+^EMA?-^)'^ . (A.3) 

i=l k=l i=l kj^l i<j k,l 

The orientifold projections are obtained by M^> = — (M^) T and Q\ 2 = sQ 2 i and 
lead to the symmetric or antisymmetric matter. Identifying Q\ 2 = Q21 leads to the 
n — 1 0{n) model. The action on M identifies for all cases = -Af , while the 
projection of the Q changes the result of the Gaussian integration starting from (A.2). 
Let us sketch the derivation of the exact loop equations for S :— Sa 2 m the following. 
The e.o.m. for the A^ eigenvalues read: 



(2) 

with an analogous expression for the X k . 

Let us define the resolvents for the two matrices as oui(z) = jj J2k' — Ho > * — 1> 2 

and choose = — ^-j, i — 1, 2 in (3.30) to obtain two Ward identities, which can be 

simplified with analogous equations to (3.31). Adding these two Ward identities one 
gets the quadratic loop equation: 

{uo^zf - + oo 2 {zf - WKzfaiz) - W^z)oo 2 {z) + h{z) + f 2 (z)) = . (A.5) 
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Here we defined 

k=l Z A k 

The derivation of the cubic loop equations is likewise very similar to the discussion 
in the orientifolded model. One uses = J^mLi (i) 1 ( 2 ) — htj i n (3.30) and sim- 
plifies the Ward identity with (3.41). Subtracting the same equation, but with indices 
(1) <-> (2) exchanged, yields the cubic loop equation: 



^{z^zY - ^{zfu^z) + W^Mz) 2 + h{z) - WKzfaiz)) - 

W^(z)Mz) 2 + h{z) - Wi(z)cj 2 (z)) + g 2 (z) - 9l (z)j = 0, (A.7) 
where 

g ^ z > ~ N 2 2^ 2^ rA « _ _ \^\ ' { ' 

k=l m=l \ A k Am )\ Z A k J 

Note that the cubic loop equation for the quiver model does not depend explicitly on 
jj. In contrast the cubic loop equation for the 0(1) model does contain a sub leading 
-^2 term and the orientifold model contains even terms at order This is consistent 
with the fact that terms come from non-orientable diagrams and these are absent of 
course in the A 2 and 0(1) model. According to the DV conjecture, terms of order are 
related to gravitational corrections. The effective pure field theoretical superpotential 
of the orbifold model can therefore be obtained by taking the effective superpotential of 
the A 2 model and simply identifying the relevant variables, e.g. = ■ However, 
in the terms related to the gravitational couplings this procedure would not give a 
correct result, since then one has to take into account the additional terms in the 
cubic loop equation! It would be interesting to investigate this issue further and to see 
if these extra terms can be related to the presence of twisted sectors in the dual string 
model. 



A. 2.1 The large A" Riemann surface for the ^-quiver 

Here we will take the large A" limit of the loop equations for the A 2 quiver. We will 
see that the Riemann surface governing the A 2 model is not hyperelliptic but rather 
is a triple cover of the z-plane, one of whose branches gives u it0 (z). The occurrence of 
Riemann surfaces with a more complicated sheet structure for multi matrix models is 
not unexpected [40, 2]. Our derivation from the exact loop equations will provide a 
proof for the curves of A 2 quiver and the orbifolds thereof. 
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The leading terms will be extracted from the loop equations (A. 5, A. 7). Using the 
variables 12 

u 1 (z)=u 1 (z)-t 1 (z) , t 1 (z) = ~(2W{(z)+W^z)) 
uj 2 (z) = -u 2 (z)+t 2 (z) , t 2 {z)= l -{Wl{z) + 2W! 2 {z)) (A.9) 
one obtains the from the large iV-loop equations the affme complex algebraic curve 

ui(z) 2 + u 1 (z)u 2 (z) + u 2 (z) 2 = p(z) 

U!(z)u 2 (z) 2 + u^zfu^z) = -q(z) (A.10) 

where 

p(z) =t\ + t x t 2 + tl-f x -f 2 

q(z) = -t\t 2 - htj + ht 2 + M - gi + g 2 . (A.ll) 
Now (A. 10) are the Viete relations for the cubic 

2 

IJ(u(z) - Ui) = u 3 - p(z)u - q(z) 

i=0 

2 

= l[(u - U(z)) + u{h + f 2 ) + 9l -g 2 - f 2 h - ht 2 = , (A.12) 

where u = —Ui — u 2 and t = —t\ — t 2 . As we have shown that /j and gi are 
polynomials in z, (A.12) is an algebraic curve, and can be viewed as a triple cover of 
the z plane [2]. The perturbations encoded in f\ and gi can be explicitly related to 
the eigenvalue densities of the two matrix model. If d\ + l,d 2 + 1 are the degrees of 
the potentials W±,W 2 , the degrees of fi,gi are di — 1 respectively. The perturbation 
parameters introduced by gi are = max(d 1 , d 2 ) and #(/) — di + d 1 for the /j. 

In particular for d x = d 2 = d relevant for the modding below we get 3d, which is is 
compatible with number of branches calculated in [36] : for equal degree d + 1 of the 
Wi it is dR + , where R + are the number of positive roots for the quiver ADE group. It 
has also been checked [36] that this number coincides with the number of independent 
S 3 in H 3 (M), where M is the local Calabi-Yau threefold after the large iV transition. 
Using the Z 2 identification \^ = —X^ on gets 

uji(z) = uj(z), uj 2 (z) = -u(-z), (A.13) 
W[{z) = W'(z), W' 2 {z) = -W'(-z) 

12 In order to keep the notation simple, we use here w, for Wj ; o- 
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etc. Since it was shown in section 3.1 that all log dependence can be absorbed into 
V(z) we can in general replace W'(z) by V'(z) and recover the large N limit of the two 
loop equations with the Z 2 orientifold identification (3.51). 
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